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Abstract 

In this paper, we study an integrable system with both quadratic and cubic nonlinearity: 
mt — bu x + \k\ [m(u 2 — u x )\ + ^k2(2mu x + m x u), m = u — u xx , where b, k\ and ki 
are arbitrary constants. This model is kind of a cubic generalization of the Camassa-Holm 
(CH) equation: m t + m x u + 2mu x = 0. The equation is shown integrable with its Lax 
pair, bi-Hamiltonian structure, and infinitely many conservation laws. In the case of b = 
0, the peaked soliton (peakon) and multi-peakon solutions are studied. In particular, the 
two-peakon dynamical system is explicitly presented and their collisions are investigated in 
details. In the case of b ^ and &2 = 0, the weak kink and kink-peakon interactional solutions 
are found. Significant difference from the CH equation is analyzed through a comparison. 
In the paper, we also study all possible smooth one-soliton solutions for the system. 
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1 Introduction 

The Camassa-Holm (CH) equation 

rat — bu x + 2mu x + m x u = 0, m = u — u xx , (1) 

was derived by Camassa and Holm [I] as a shallow water wave model. In recent years, this 
equation has attracted much attention in soliton theory. In the literature, this equation was 
implicitly implied in the paper of Fuchssteiner and Fokas on hereditary symmetries as a very 
special case [2J. However, since the work of Camassa and Holm pQ, various studies on this equa- 
tion have been remarkably developed. The CH equation possesses many important integrablc 
properties. For instance, it admits Lax representation, bi-Hamiltonian structures, multi-soliton 
solutions, and algebro-geometric solutions [3]-[S]. Also, it is integrable by the inverse scattering 
transformation The most remarkable feature of the CH equation ([1]) is to admit peaked 

soliton (peakon) solutions in the case of b = 0. A peakon is a weak solution in some Sobolev space 
with corner at its crest. The stability and interaction of peakons were discussed in several refer- 
ences [IS] -[IE]- I* 1 addition to the CH equation being an integrable model with peakon solutions, 
other integrable peakon models have been found. Those models include the Degasperis-Procesi 
equation [H]-[22] and the cubic nonlinear peakon equations [23j-[26j. 

The present paper focuses on the following equation with both quadratic and cubic nonlin- 
ear ity: 

m t = bu x + ^k 1 [m(u 2 -ul)] x + ^k 2 (2mu x +m x u), m = u-u xx , (2) 

where b, k\ and k 2 are three arbitrary constants. It is clear that equation ([2]) is reduced to the 
CH equation ([I]) when we take k\ = 0, k 2 = —2. For k\ = —2, k 2 = 0, equation ([2]) is exactly 
the cubic nonlinear equation: 

m t -bu x +[m(u 2 -ul)] x =0, m = u-u xx , (3) 

which was derived independently by Fokas [27] . by Fuchssteiner [28], by Olver and Rosenau [I], 
and by Qiao \23\ [29] where the equation was derived from the two-dimension Euler system, and 
Lax pair, the M/W-shape solitons and cuspon solutions were presented. Recently, Gui, Liu, 
Olver and Qu worked out the peakon stability for the cubic nonlinear equation ([3]) in the case 
of 6 = [30]. 

Equation ([2]) is actually a linear combination of CH equation (pQ) and cubic nonlinear equa- 
tion ([3]). Therefore, we may view equation ([2]) as a generalization of the CH equation, or simply 
call equation ([2]) a generalized CH equation. This structure is very similar to the one in deal- 
ing with the Gardner equation, known as a linear combination of KdV and mKdV equations, 
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which has important applications in various areas of physics [3H 32J. In fact, the structure of 
Gardner equation is our mainly starting point to study equation ([2]). We also notice that by 
some appropriate rescaling, equation ([2]) might implicitly be implied in the paper of Fokas and 
Fuchssteiner in the context of hereditary symmetries [27, [28] , 

In this paper, we first present the Lax representation, bi-Hamiltonian structure and infinitely 
many conservation laws for equation ([2]). This indicates that equation ([2]) is completely inte- 
grable. Then we show that this equation possesses the single peakon of traveling wave type 
as well as multi-peakon solutions in the case of b = 0. In particular, for the one-peakon case, 
we find that the complex peakon (i.e. pekaon with complex coefficient) appears by properly 
choosing the parameters k\ and hi. For the two-peakon case, the peakon dynamical system is 
explicitly presented and their collisions are discussed in details. Significant difference between 
the CH equation ([T]) and the generalized CH equation ([5]) (for the case of b = 0) is analyzed 
through a comparison. 

Another purpose of this paper is to show that equation ^ with ki = (namely, cubic 
nonlinear equation ([3])) allows the weak kink solution in the case of b ^ 0. Different from the 
multi-peakon solutions in the form of linear superpositions of the single-peakon, equation ([2]) 
with ki = and b ^ does not allow the multi-kink solution in the form of the superpositions 
of single-kink. However, we find that equation ([2D with ki = and b ^ allows the solutions 
in the form of the superpositions of single-kink and multi-peakon. In particular, the weak kink 
and kink-peakon interactional solutions are shown and plotted. Within our knowledge, this is 
probably the first time discussing the weak kink and kink-peakon interactional solutions. In the 
paper, we also study all possible smooth one-soliton solutions. 



2 Lax pair, bi-Hamiltonian structure and conservation laws 

Let us consider the following pair of linear spectral problems 




U 



V 



u 



V 




(4) 
(5) 
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where 

■-y/l - X 2 b, 



a 



-2„. , «. / .2 ..2. 1 



(6) 



A =A~ a + —k\(u - u x ) + -k 2 (au - u^), 
B = — A _1 (u — atfj) — t^ 771 [^i( u2 ~~ n ir) + ^2^] > 
C =A _1 [fei(n + au x ) + fo] + ^ \k\m(u 2 — u 2 ) + k\k 2 {mu + u 2 — u 2 ) + fcfu] • 
One can easily see that the compatibility condition of ^ and j5]) is 

U t -V x + [U,V] = 0. (7) 

Substituting the expressions of U and V into (JT)), we find that ([7]) is nothing but equation ([2]), 
namely, (jU) and §5§ are the Lax pair of 

A direct computation shows that equation (|2|) has the following bi-Hamiltonian structure 

SHi m 2 

m t = J — — = K— — , 8 
dm dm 

where 

r+00 

J = k\dmd~ x md H — k2(dm + md) + bd , H± = - (u 2 + u 2 )dx, (9) 

2 2 y_ 00 

if = a -9 3 , H 2 = - (hu 4 + 2k lU 2 u 2 x - -k lUx + 2k 2 u 3 + 2k 2 uu 2 x + 4bu 2 )dx.(W) 

Let us now construct conservation laws of equation ([2]). Let a; = then w satisfies the 
following Riccati equation 

w x = —-\(kxm + fe) + oluj - ^Xmoj 2 . (11) 
Based on (j4|) and ([5]), we obtain 

{ln<j)i) x = ~ + -\moj, Qjx<f> 1 ) t = ~A--Bu, (12) 
which yields the conservation law of equation ([2]) : 

Pt = F x , (13) 

where 
p = mw, 



F = -m \k\{u 2 — u 2 ) + k 2 u\ lj — - [afci(u 2 — u 2 ) + k 2 (au — u x )] A -1 + (ti — oiu x )oj\~ 2 . 



(14) 
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Usually p and F are called a conserved density and an associated flux, respectively. Next, we 
derive the explicit forms of conservation densities in the case of b = 0. In this case, equation 
(fTTl) becomes 



a:,. — -\{kim + k 2 ) + ui — -Amu; 2 



2 v w 2 

We expand co in terms of negative powers of A as below: 



(15) 



j=0 

Substituting (fTB"j) into (fT5j) and equating the coefficients of powers of A, we arrive at 

,2 



(16) 



UJ j+1 



1 / — ; 5 ; 2k\m l + A; 2 (2m + m a 

— y , -k 1 m z -k 2 m, ui = 3—— „— 



m 



1 



ujj - uj jtX - -m 



E 



<jji<jjk 



(17) 



, J>1- 



Inserting (|16p and (|17p into (|14p , we finally get the following infinitely many conserved densities 
and the associated fluxes 



Po 
Pi 

Pj+i 



\J — k\m 2 — k~2m, Fq = — [ki(u 2 — u\) + k 2 v] \/—kim 2 — k 2 m, 

2kim 2 + k2(2m + m x ) 1 r 2 2 -, 1 r 2 2 i 



1 / 1 

r' 1 u * 2" 



i+fe=j"+l, i,fc>l 



(18) 



Fj+i = (u - «a,)a;j_i + - [fci(u 2 - u 2 ) + k 2 u] p j+1 , j > 1. 

We can also consider the expansions of oj in the positive powers of A: 



3=0 



Substituting (fT9j) into equation (fT5j) and comparing powers of A lead to 



W2j =0, j > 
1 
2 



(19) 



= ^[fa-fa + n z) + fa], u 2 j+i,x ~ ^2j+i = ~m ^2 Via*., j > 1. 



(20) 



Equation (|20j) shows that one may solve a first-order differential equation to obtain nontrivial 
W23+I1 J > 1 in this sequence. For brevity, we omit the result since this may involve in nonlocal 
expressions in u. 
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3 Peakon solutions in the case of b = 

Applying the operator (1 — d x )~ x to equation ([2]), we obtain 

ut =\ki(3u 2 u x - ul) + \k 2 uu x + \ki{l - dl)" 1 ^ 

1 / 2 1 > < 21 > 

+ - dly 1 ( 26u + fci(uu* + -u 3 ) + fc 2 (u 2 + -u x ) 

Taking the convolution with the Green's function G(x) = ^ exp(— \x\) for the Helmholtz operator 
(1 — d x ), equation (|2T|) can be rewritten as 

u t =^ki(3u 2 u x - ul) + ^k 2 uu x + ^kiG(x) * u 3 x 

In this section, we derive the single and multi-peakon solutions of equation ([2]) with 6 = 0. 
3.1 Single-peakon solutions 

Let us suppose the single-peakon solution of equation (|2|) with b = in the form of 

u = Ce-I*-*!, (23) 



where C is to be determined. The derivatives of expression (123D do not exist at x = ct, thus 
(|23p can not satisfy equation (|2|) with 6 = in the classical sense. However, in the weak sense, 
we can write out the expressions of u x , ut and m with help of distribution: 

u x = -Csgn(x-ct)e- lx - ctl , u t = cCsgn(x - cf)e -|a! - rf| , m = 2C5(x - ct). (24) 



Substituting (|23|) and (|24jl into the weak form (|22|) with 6 = 0, we are able to find that C should 
satisfy 

l kl C 2 + h 2 C + c = 0. (25) 

For k\ = 0, ^ / 0, we obtain C = In particular, when k 2 = —2, (|23l) is exactly the 

single-peakon solution u = ce~ \ x ~ ct \ of the CH equation (JTJ) with 6 = 0. For k\ = —2 and k 2 = 0, 
we recover the single-peakon solution u = ±^J^e~\ x ~ ct \ of the cubic nonlinear CH equation (|3|) 
with 6 = [30]. In general, for k\ ^ 0, we may obtain 

3 (V3k 2 ± y/3k% - I6hc 



C = * = '-. (26) 

4^3/ti 
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If 3A£ - 16ifeic > 0, then C is a real number. For example, we choose k\ = k 2 = —2 and c > 0, 
then the corresponding single-peakon solution is u = ~ 3=l=v ^ 9+2 — e~l x ~ d L If 3/c| — 16&ic < 0, then 
C is a complex number. This means the peakon solution with complex coefficient is obtained. 
For example, we set A>2 = —2 and k% = -, then we have u = ( 3± ^ t ) c e ~l :r ~ ct l. 

3.2 Two-peakon solutions and their dynamics 

Let us assume the two-peakon solution to equation ^ with 6 = has the following form 

u = Pl (t) e -I*-»WI +p 2 (*)e" |a, " aa(t)l . (27) 

The expression of u has two peaks (or troughs) at positions x = qi(t) and x = q 2 (t). u x , Ut and 
m have the following weak forms: 

u x = — pisgn(x — qi)e~\ x ~ qi \ — p 2 sgn(x — q2)e~^ x ~ q2 K 

u t = Vi,te~ ]{X ~ qi]{ +P2,te" |x '" 921 +piqi it sgn(x - gi)e" |x '"' ?l1 + P2q2,tsgn(x - q 2 )e~ ][X ~ q2 \ (28) 
m = 2pi5(x - qx) + 2p 2 5(x - q 2 ). 



Substituting (|27p and f|28|) into (|22p with 6 = 0, we may obtain the following peakon dynamical 
system 



Pi,t = -\k 2 p\p 2 sgn{q\ - q 2 )e \ qi 92 l, 

(29) 



P2,t = -\k 2 piP2sgn(q2 - q\)e 192 



9i,t = -k lPl p 2 e^-^\ - \kxpl - \k 2 (pi +p 2 e-^' q ^) , 
q2,t = -kipip 2 e-\^- q2 \ - \kxpl - \k 2 (p 2 + Pie^ 1 ^) . 

For k\ = 0, ^2 = — 2, ([2]) recovers the CH equation and (|2"9"j) is reduced to the two-peakon 
dynamics of CH equation with 6 = 0, which Camassa and Holm studied in PQ. In particular, 
the peakon-antipeakon was taken as below 

u = Pl (t) ( e -l^(*)l - e H*+9i(*)l) , 

where 

Pl (t) = ccoth(ct) = -p 2 {t), qi (t) = - In — + 2 e _ rf = -q 2 {t). (30) 

Camassa and Holm [T] also pointed out that the collision of the peakon-antipeakon occurs at 
the moment t = 0. At the moment of collision, the amplitudes p\ and p2 become infinite, but 
the solution u converges to zero [IJ [33] . 

For k2 = and 6 = 0, ([2]) is exactly the cubic nonlinear equation [U [231 G3 EE] 

m t = -kx [m(u 2 - u 2 x )] x , m = u- u xx , (31) 
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and the peakon dynamical system (|29p is simplified to 

Pl,t = 0, P2,t = 0, 
< ?i,t = -kw&e-fo-rt ~ \k x p{ , (32) 
k 92,t = -fcipi^e" 191 " 921 - gfeipi- 

From the first two equations of (f32j) . we obtain = ci, P2(i) = C2, where ci, C2 are two 

arbitrary constants. For the case of c\ = c|, the solutions of and (72 (i) in (f3"2"j) are given by 

gi (i) = -fci(c lC2 e-l^-^l + ±c?)t + Ax, 
g 2 (t) = -fc 1 ( Cl c 2 e-l Al -^l + ±c?)t + A 2) 

where Ai, A2 are arbitrary integral constants. If A\ = A2, the two-peakon solution is reduced 
to zero (when C2 = —c\) or a single-peakon (when C2 = c\). li A\ ^ A2, the two-peakon can not 
collide since q\(t) 7^ cti(t) for every i G (—00, +00). Especially, in the peakon-antipeakon case 
(ci = — C2), the collision can't occur for Ai 7^ A2. This is very different from the case of the CH 
equation [I]. For the case of cf / c|, we obtain the following solution: 

' Pl(*) = ci, P2(*) = c 2 , 

< gift) = ^(0 | fc ?(|i C |)| fe-li* 1 ^-^' - l) - \k^t, (34) 
^ Q2(t) = sgn(t)^m^ (e-ll^^-D*! - l) - ifci^t. 

The two-peakon collision occurs at the moment t = 0, since <?i(0) = 92(0) = 0. Without loss 
of generality, let us suppose cf < cf. From the formula (f3"3| . we know that for i < 0, the tall 
and fast peakon (with the amplitude |c2 1 and peak position g 2 ) chases after the short and slow 
peakon (with the amplitude |ci| and peak position q\). At the moment of t = 0, the two-peakon 
collides and overlaps into the peakon: u(x,0) = (ci + C2)e~' a; '. After the collision (t > 0), the 
two-peakon departs, and the tall and fast peakon surpasses the short and slow one. See Figures 
1 and 2 for the developments of this kind of two-peakon. 

Remark 1. Our results show that the collision of two-peakon of equation (|31[) is very 
different from the case of CH equation. In the peakon-antipeakon case, the collision can occur 
for the CH equation, but could not happen for equation (|3ip . The collision of two-peakon of 
equation ([3T|) occurs in the "chase" case as shown above (see Figures 1 and 2). 

Next, we want to derive the solutions of (|29p for the general case k\ 7^ and k,2 7^ 0. To do so, 
let us introduce the transformations p{t) = p\{t) +P2(t), q(t) = qi(t) + g 2 (i), P(t) = p\(t) — P2(t), 
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Q{t) = Qi(t) ~ Q2(t)- Then (|29|) is transformed to 
Pt = 0, 

qt = -^h(p 2 - P 2 )e-IQI - lh(p 2 + P 2 ) - \k 2 p{l + e-WI), 



^ = -|A; 2 (p 2 -P 2 ) S5 n(g)e-l«l, 
_ Q t = ±& 2 P(e HQI - 1) - ifcipP. 

From the first equation of (|35j) . we know that p(i) = Ai, where Ai is an arbitrary constant. Let 
r = 1 + I^Aij then we have the following results. 

Case 1. If < T < 1, then (|35|) admits the following solution 



P(t) = ±02 i±4<!, 0{t) = ±ln 



1-A 3 e m a 2 (l + A 3 e Bt ) 2 - A 2 {1 - A 3 e Bt ) 2 

, l^3e B '-^±f|| 2fc ia 2 (3r-l) 1, 4 1 „ 

o(t) = - In— — ~ - - AfcAi + -fci(i4? + aS) i + A 

|A 3 e Bi -^f|| 3P(A 3 e Si -l) 2 L 3 n 1 2J1 

where a 2 > |^i|, A 3 > 0, B = — ic^for, and P is an arbitrary constant. 
Case 2. If T > 1, then ([35]) has the following solution 



F(() = ±02 i^, 0{t) = ±ln 



These yield 



3(e 2 -t T ) +ln(e24 + 1) 4 1 



= coth(t), qi (t) = ol ^ t Ts + ln(e 2t + 1) - £t- In 2 

5 



Pa(t) = - coth(t), g 2 (t) = — i _ - - ln(e 2t + 1) + |t + In 2. 



(36) 



(37) 



l + i 3 e Bi ' a 2 (l -i 3 e B ') 2 - A 2 (l + A 3 e Bt ) 2 ' 

i l^3e B ' + ^f|| 2/fc ia 2 (3r-l) l r , , 1 ,, 2 2 ,. „ 

a(t) = - In—, ^ 2 H ^ feMi + -fci(A 2 + a?) i + P, 

l^e Bi + fe| 3P04 3 e^ + l) 2 L21 3 Ul 2)1 

where < a 2 < ^=-, ^3 > 0, P = — ^a 2 /c 2 r, and P is an arbitrary constant. 
Let us now give special two-peakon solutions based on Case 1 and Case 2. 
For Case 1, let Ai = 0, a 2 = 2, A 3 = 1, k\ = fc 2 = —2, and P = 0, then we have 

16 4 4e 2 * 

p(t) = 0, P(i) = 2 coth(t), g ( t ) = __ + -t, Q(t) = - ]n __ . (38) 



(39) 



Therefore, we obtain the following peakon-antipeakon solution 

u(x, t) = coth(t) (e-l^-^WI - e-1 x -»Wl) , (40) 
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where q\ (t) and 92 (i) are shown in (|39p . From (|38p . one can easily know that the collision occurs 
at the moment t = since Q(0) = 0. From (I39p . we may compute 

limpi(i) = — limp2(*) = 00, ]imqi(t) = lim (72(f) = 00. (41) 
t— >o t->o t— >o t— >o 

But from ([40 p . we may infer that 

lim u(x, t) = 0, for every x £ R, (42) 
t-fO 

which indicates that the peakon and the antipeakon vanish when they overlap. Guided by the 
above results, we may describe the dynamics of peakon-antipeakon solution (|40p as follows. For 
t < 0, the peak is at q% (t) and the trough is at q±(t). The peak and the trough approach each 
other as t goes to 0. At the moment of t = 0, the peakon and the antipeakon collide and vanish. 
After their collision (t > 0), they depart and redevelop with the trough at qi (t) and the peak at 
q\{t). Figure 3 shows the peakon-antipeakon dynamics. 

Remark 2. The amplitudes p±(t), P2(t) in formula f|39|) are the same as those of the CH 
equation p], but the peak positions qi(t), q2(t) are different (see formulas ([30]) and ([39|) ). In 
the CH equation, only pi(t) and pzit) become infinite at the instant of collision [HE]- I n the 
new equation ([2]) with 6 = 0, not only pi(t), P2(t) but also qi(t) and q2 (t) become infinite at the 
instant of collision (see formula (|4ip ). However, in both case, the peakon-antipeakon vanishes 
when the overlap occurs (see formula (|42p ). 



For Case 2, choosing ki = &2 = —2, A\ = 3, A3 = 1, T = 3, 02 = \/3, B = 3\/3 forces the 
first two equations in ([37|) to 

_ e ?j\f?>t g g 3\/3t 

p v = ^®T^#m> Q ^ = S9n ^ ln e6 ^ + 4e 3^ + i - (43) 

Substituting ([4"3"|) into the second equation of ([33]) . we may obtain the expression of g(t) and 
finally arrive at the solution 

, . 3 y/3 , , ,Sy/3t. . . 3 \/3 , , ,3\/3i N 
= 2 ~ -y s 5 f "(i)tanh(^— ), p 2 {t) = - + — s£in(t) tanh(— ^— ), 

1 6e 3v ^* 1 ^ e 3v/ 3* + 2 + \/3 16\/3 7 

^ 1 6e 3v ^ 1 e 3v/3 ' + 2 + \/3 16^3 7, 
92 1 = sqn(t) in = = In = — = — t. 

2 y eSv^i + 4e3>/at + 1 2 e 3^3t + 2-^3 9(e 3 ^ 3t + 1) 2 

The formulas (|43p show that the collision happens at the moment t = since Q(0) = 0. 
From (|44p. we know that the two-peakon collides and overlaps into the peakon u(x, 0) = 
3 e -k+2 ln ( 2 +v / 3)H g- 1 a ^ ^-^g momen t ^ = 0. After the collision, the two-peakon departs and 
redevelops. See Figure 4 for the profile of the two-peakon dynamics. 
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3.3 N-peakon dynamical systems 

In general, we suppose that iV-peakon has the following form 

N 

u(x,t) = Y,Pj(t)e~ lx ~ qm . (45) 
i=l 

Similar to the cases of one-peakon and two-peakon, but with a lengthy calculation, we are able 
to obtain the following 7V-peakon dynamical system 

1 N 

Pj,t = -^Pj^Pksgniqj 

k=l 
1 - 
k=l 

Different from the iV-peakon system of the CH equation, the above system can not directly be 
rewritten in the standard form of a canonical Hamiltonian system. It is very interesting to see 
whether the above system is an integrable Hamiltonian system under a Poisson structure. We 
will investigate this in the near future. 

4 Weak kink solution in the case of b ^ 

We have already shown that equation ([2]) admits peakon solutions in the case of b = 0. It is 
natural to ask what kind of solution one may obtain for the case of b 7^ 0. Here we will reveal 
that equation ([2]) with k<i = and 6/0 (i.e. cubic nonlinear equation (J3|) with b / 0) possesses 
the weak kink and kink-peakon interacted solutions. 

4.1 Single weak kink solution 

We seek the solution of equation ([2]) in the form of 

u = Csgn(x - ct) (e _ - l) , (47) 
where the constant C is to be determined. If C 7^ 0, (|47p actually is a kink solution since 

lim u = — lim u = —C. (48) 

x— >-+oo x— >— 00 

One may easily check that the first order partial derivatives of (|47p read 

u x = -Ce-\ x ~ ct \ ut = cCe-\ x ~ ct \. (49) 



Qk)e 



Mj-ik\ 



1 ,/l 2 f , 

V i,k=l 



sgn{qj - qi)sgn(qj - q k ))e 
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u(x,t) 

3- 

2.5 1 , 
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Figure 1: The two-peakon solution determined by 
([M]) with ci = 1, c 2 = 2, fci = -2. Black line: £ = 
—4; Red line: £ = — 1; Brown line: £ = (collision); 
Blue line: £ = 1; Green line: £ = 4. 

u(x,t) 




Figure 3: The peakon-antipeakon solution (|40|) . 
Pink: peakon (and antipeakon) with the peak (and 
trough) position qi\ Green: antipeakon (and peakon) 
with the trough (and peak) position q\. 




u(x,t) 




Figure 2: The two-peakon solution determined by 
with C = -1, c 2 = 2, fci = -2. Black line: £ = 
—6; Red line: £ = —2; Brown line: £ = (collision); 
Blue line: £ = 2; Green line: £ = 6. 

u(x,t) 
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Figure 4: The solution u(x, t) with pi, p2 and qi, 
(72 shown in (|4^1) . Black line: £ = —6; Red line: 
£ = —0.8; Brown line: £ = (collision); Blue line: 
£ = 2; Green line: £ = 4. 

u(x,t| 




Figure 5: The kink solution ([54]) for £ = 0. 



Figure 6: The kink-peakon interacted solution. 
Black line: £ = 2; Blue line: £ = 1; Green line: £ = 0. 
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The second order partial derivatives of ()47|) do not exist at x = ct. Therefore, like the case of 
peakon solutions, the kink solution in the form of (I47p should also be understood in the weak 
sense. 

Substituting (|3T|) and (|4"9j) into the weak form (j22|) . we arrive at 

cCe -\x-ct\ = \[k 2 C 2 (x - ct) - bC]e-\ x ~ ct \ - ^[hC 3 + bC](x - ct)sgn(x - ct^*-^ . (50) 



If two sides of equation (|50j) match, we should require 

' k 2 = 0, 
-\bC-cC = Q, 
hC 3 + bC = 0, 



(51) 



which leads to 



k 2 = 0, 

C = 



\b, 



(52) 



Remark 3. Formula (|52|) shows that equation ([5]) possesses the weak kink solution in the 
form (|47p for the case of k 2 = and 6 7^ 0. In formula ()52|) . k 2 = implies that the CH equation 
(prj can not allow a weak kink solution in the form of (|47p . c = —^b means the kink wave speed 
is exactly —\b. This is very different from the single-peakon solution (the wave speed in the 
single-peakon (I23D is an arbitrary constant c). 



m = u — u-, 



In particular, we take k 2 = 0, k\ = — 6 = 2, then equation ([2]) is cast to 
m t + 2u x - [m(u 2 - u 2 x )] x = 0, 
and the corresponding weak kink solution is 

u = sgn(x - t) ^e"'^*' - 1 
See Figure 5 for the profile of this kink wave solution. 



(53) 



(54) 



4.2 Kink-peakon interacted solutions 

Let us first point out that equation ([2]) with b 7^ does not allow the two-kink solution in the 
form of the superposition of two single-kink solutions: 



u = Pl {t)sgn{x - qi (i)) (V |x - 9lWI - l) + P2 (t)sgn{x - q 2 [t)) (V |x -" 2(t)l 



(55) 
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In fact, substituting (|55p into the weak form (|22p of equation ([2]), we find that the solution 
assumed in the form (I55p is reduced to zero or single-kink solution (|47H . 

Instead of considering the solution in the form (|55p . let us make the following ansatz of 
solution 



u = Pl (t)sgn(x - qi (t)) (Vl*-9iWI - l) + p 2 (t)e^ x - q2 ^ , 



(56) 



which actually describes a new phenomena of kink-peakon interacted dynamics in soliton 
theory. Apparently, with the help of distribution, u x and ut can be calculated below: 

u x = -pie~\ x ~ qi \ -p 2 sgn(x - q 2 )e~^ x ~ q2 \ 

u t =p u sgn{x- qi {t)) ( e -l*-9i(t)l - l) +p 2j , e -l^l + pm>te -\x- q i\ + p2 q 2tSgn ( x _ g 2 ) e -l-^l. 
Substituting (|56p and (|57p into (|22p with = 0, we arrive at 



(58) 



P2,t = kiPiP2sgn(q 2 - qi)e~\ qi ~ q2 \ , 
qi,t = ~\b ~ kiPiP2sgn(q 2 - qi)e-\ qi ~ q2 \ , 
I, <3V = -\k\p\ - \k\p\ + h [pi - pip 2 sgn(q 2 - qi)) e~ l<?1 ~ 921 + k 1 sgn(q 2 - qi)pip 2 . 

Let us choose k\ = — b = 2, then pi = dbl. Without loss of generality, taking p\ = 1 forces (f58|) 
to 

P2,t = 2p2sgn(q 2 - qi)e~\ qi ~ q2 \ , 

q ltt = 1 - 2p 2 sgn{q 2 - g x ) e -l«n-«al , (59) 
. ?2,t = - 1 + 2 (1 - p 2 sgn(q 2 - qi )) e -l«-*»l + 2p 2 sgn{q 2 - qi ). 

To solve the above system, let us make an assumption qi < q 2 . Therefore, integrating equation 
(f59l) . we may obtain 

qi = t-p 2 + At, 

q2 = t -p 2 -\n(\pl- ipa + i + ^ + Ai, (60) 

I P2,< = |P2-P2 + 2 P2 + -42, 

where j4i and y4_2 are two arbitrary constants. Letting A 2 = 0, then we may solve the third 
equation of (|60p for p 2 , which has the following implicit form: 

1 . o 9 _ 3V7 



In \p 2 



In 



;P2 + 9 



4»9 — 9 

+ arctan( — ) = 2t + A 3 . 
3v 7 



(61) 



See Figure 6 for the profile of the kink-peakon interactional solution with A\ = A 2 = A3 = 0. 
Here, we do not have more further explicit results about the collisions between kinks and peakons 
due to the implicit form of the kink-peakon solution. 



15 



In general, we may make the following ansatz of the solution for equation (|2j) with k 2 = 
and 6/0: 

N 

u = Po (t)sgn(x - q (t)) ^ x ~ qo ^ - lj + ^p j {t)e~\ x -' 1 ^ , (62) 

i=i 

which can be viewed as the interaction of single-kink and iV-peakon solution. Through a very 
lengthy calculation, we are able to obtain the following interactional dynamical system of single- 
kink and iV-peakon: 



i n 

-kipl + hpo^pisgniqo - qi)e l<?0 9l1 



st 2 

i=l 



+ 2 fc l ^2 PiPksgn(qi - qk)(sgn(q k - q ) - sgn(qi - q ))e ^ 
i,k=l 

N 

Pj,t = kipoPjsgn(qj - go)e~' 9o ~ 9il + hpoPj^Pisgnfaj - qi)sgn(qj - q )e~^~ qi \ 



(63) 



i=i 

11 1 N 

qj,t = gfci^ 2 - -k lP 2 (l - 2e^ qo ~ q ^) - -h PtPfc(l " s<M<7j - %)^n( gj - g fc )) e -k>-*H*-**l 

i,k=l 

N 

- hpo^Pi (sgn( qj - q )(e^ q °- q ^ - l) e -l*-*l - sgnfa - gi )e-l»-*H»-fcl 
i=i 

The above system is not presented in the canonical Hamiltonian system. We still do not know 
whether this system is integrable under a Poisson structure. 



5 Smooth soliton solutions 

In this section, we study all possible single smooth soliton solutions to equation ([2]). To do so, 
let u(x,t) = 4>{x — ct) = </>(£), where c is the wave speed. Substituting it into equation (|2|) and 
integrating it, we have 

f (c+^W+^ 2 -^iW')j = ((& + c)</>+^M> 2 + ^i/) " Q*2 + ^M) (0') 2 ,(64) 

where "/" stands for the derivative with respect to £. Equation (|64p can be rewritten as the 
following two-dimensional system 

d4 = dy_ - {\k 2 + \k^) y 2 + + c) + \k 2 $ + \k x <\> 2 ) 

di v " di c + \k 2 ct> + ±m 2 - \k iy 2 ' [ ' 
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which has the following first integral 

H(<f>, y) = \y 2 (c + \k 2 4> + - h iy * - ((6 + c)0 2 + + = h. (66) 

Without loss of generality, we assume that the wave speed c is a positive constant. Then, 
system ([65]) is a three-parameter planar dynamical system with the triple tuple (b, fci,/^)- We 
only pay attention to the bounded solutions of ([65]) . 

We notice that for k\ ^ 0, the right hand of the second equation in (|65p is not continuous 
on the hyperbola c + ^24> + \ki<f) 2 — \k\y 2 = 0, i.e., [6 + Jrj^J — y 2 = ki . In other words, 
<f>g is not well-defined on this curve in the phase plane (<j>, y). System (|65p is called the second 
kind of singular travelling wave system [M] . 

Let us assume c > 0, k\ ^ 0. Imposing the transformation cZ£ = (c+ 5^2^ + \k\<\) 2 — \k\y 2 )dC, 
for c + \ki§ + \k\(\> 2 — \k\y 2 / on system ([6"5]) leads to the following cubic system 



f = y( c+ i^+iA: 1 2 -ifc 1 y 2 ), 
I = - + y 2 + H(b + c) + |fc 2 + ±M 2 ), 

which can be rewritten as 

§L = -vy 2 + ^ + A^ + B = -w 2 + ^-^) ^ + ^ + ^0d4 



(67) 



(68) 



where <f> = <p - fc, r, = R = Mz*£, A = ^(ffi"* , and 5 = MMzfW£» , Apparently, 
(1681) has the following first integral: 

flifoi/) = - y 2 ) 2 - ~ V - -B V = h. (69) 

Thus, we have 



.2 



(y? 2 - R) ± V(^ 2 - 4/i) - 4Bp - 2p + R)ip 2 . (70) 



To investigate the exact soliton solutions of equation ([2D, let us consider system ([68]) with 
B = 0, namely homogeneous system. 

1. M-shape and W-shape soliton solutions 

Assume k\ < 0, k 2 = Ak\{b + c) > 0, b + c < 0, and let wi = ^tjV^ 2 — c 2 . Then, system 
(]68p is able to be solved with the following smooth M-shape and W-shape soliton solution for 
equation ([2]): 

. (b-c)J%g£ cosh( X ) 
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When — oo < b < —3c, the functions defined by (|7ip take on the graph of smooth M-shape 
and W-shape profiles shown in Fig. 7 (7-1). When —3c < b < — c, the functions defined by (I71D 
give smooth solitons shown in Fig. 7 (7-2). 




(7-1) 99(0- (7-2) </;(£). 

Fig. 7 The profiles of waves with respect to £ 

2. Smooth kink and anti-kink wave solutions 

Suppose k\ > 0, h\ = 4fci(6 + c) > 0, 6 + c > 0, 6 < 0, c > 0, and let 02 = \/jS:> w 2 
^-a/c(6 + c). Then, system dggj may be solved with the following kink and anti-kink wave 
solutions for equation ([2]): 



fix) 



y/\b\(b+c) sinh(x) 



(^+ v / |b|cosh( X )' 



2c X ~ glr ln TTO? + V 7 ^^) In (-+D("l+-+^v / ^cos hx ) 



(72) 



(to— l)(a|— w+«2 V'I+ckI coshx) ' 

where u; = a2sinh(x). 

However in the case of &i > 0, A;| = 4fci(6 + c) > 0, 6 = 0, we can solve the system (|68|) with 
the following kink and anti-kink wave solutions for system ((21) in a simpler form: 



v?(x) 

e(x): 



1+e-x ' 

±|(X + In coshx)- 



(73) 



3. Smooth soliton solutions 

Suppose ki > 0, fc| = 4fci(6 + c) > 0, < b < c, and let 03 



b -, and u 3 



Then, we can solve the system (|68p with the following soliton solutions for equation ([2]) 



fix) 
Six)- 



c-b 



' v / c+v / 5cosh(x) ' 



2 <* ~ ^ In ^fgg + V / ^b) In (^X-i+^v^oshx) 



(74) 



(to— l)(a|— w+a3v/l+a| coshx) ' 
ct3sinh(x). The functions defined by (J7JJ) present two soliton solutions shown in 



where w 
Fig.8 (8-1). 

4. Two-crest soliton solutions 
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Suppose k\ > 0, &2 = 0, < b < c or b > c > 0. Then, we can solve the system ([68]) with the 
following two-crest soliton solutions for equation ([2j) 




(8-1) <p(0- (8-2) <p(£). 

Fig. 8 The profiles of soliton waves with respect to £ 

Remark 4. If B ^ in system (|68p . namely inhomogeneous system, then we may still 
obtain some exact smooth one-soliton solutions to equation ([2]) |35j . 

6 Conclusions and discussions 

In this paper, we have presented the Lax representation, bi-Hamiltonian structure and infinitely 
many conservation laws for equation ([2]). The peakon solutions for this equation are derived 
in the case of b = 0. For the one-peakon solution, we found that the complex peakon can be 
obtained through properly choosing the parameters k\ and ki in equation ([2]). For the two- 
peakon solution, its peakon dynamical system was explicitly solved and their collisions were 
shown in Figures 1, 2, 3, 4. Moreover, we found that equation ([2]) with ki = (namely, cubic 
nonlinear equation ([3])) possesses weak kink and kink-peakon interactional solutions in the case 
of b ^ 0. Compared with the CH equation, equation ([2]) has some different features, such 
as the dynamics of the two-peakon solutions and the existence of the weak kink and kink- 
peakon interactional solutions. These differences are mainly caused by both quadratic and 
cubic nonlinearity in equation ([2]). In the paper, we also study all possible smooth one-soliton 
solutions. Other topics, such as cuspons, Darboux transforms, peakon stability, and algebra- 
geometric solutions, remain to be developed. 
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